Hadronic centrality dependence in nuclear collisions 



Sonja Kabana 

Laboratory for High Energy Physics, University of Bern, Sidlerstrasse 5, 3012 
Bern, Switzerland, e-mail: sonja. kabana@cern. eh 

Abstract. The kaon number density in nucleus+nucleus and p+p reactions 
is investigated for the first time as a function of the initial energy density e 
and is found to exhibit a discontinuity around e=1.3 GeV/fm 3 . This suggests 
a higher degree of chemical equilibrium for e > 1.3 GeV/fm 3 . It can also be 
interpreted as reflection of the same discontinuity, appearing in the chemical 
freeze out temperature (T) as a function of e. The N "" 1 dependence of (u,d,s) 
hadrons, whith N the number of participating nucleons, also indicates a high 
degree of chemical equilibrium and T saturation, reached at e >1.3 GeV/fm 3 . 
Assuming that the intermediate mass region (IMR) dimuon enhancement seen by 
NA50 is due to open charm (DD), the following observation can be made: a) 
Charm is not equilibrated, b) J/ty/DD suppression -unlike J/fy/DY- appears 
also in S+A collisions, above e ~1 GeV/fm 3 . c) Both charm and strangeness 
show a discontinuity near the same e. d) J/W could be formed mainly through cc 
coalescence, e) The enhancement factors of hadrons with u,d,s,c quarks may be 
connected in a simple way to the mass gain of these particles if they are produced 
out of a quark gluon plasma (QGP). We discuss these results as possible evidence 
for the QCD phase transition occuring near e ~1.3 GeV/fm 3 . 



1. Introduction 

The quark-gluon plasma phase transition predicted by QCD |jj may occur and 
manifest itself in ultrarelativistic nuclear collisions through discontinuities in the 
initial energy density (e,) dependence of relevant observables. A major example of 
a discontinuity is seen in the J/^/DY [|| discussed e.g. in ||, ||. We investigate here 
for the first time the dependence of strangeness production, in particular of kaons, on 
the initial energy density e^. The degree of equilibrium achieved in nuclear collisions 
has been intensively studied comparing hadron ratios and densities to models (see e.g. 
1@' H 0)' ^ e mves tigate here if chemical equilibrium is achieved, examining an 
other aspect of equilibrium states, namely the volume (V) independence of hadron 
densities (p). 

2. Results and discussion 

The kaon density (/9r-=(K per collision)/ V) at the thermal freeze out in nuclear 
reactions, investigated as a function of the initial energy density t\ (figure l,(a)) (see 
Q for calculation details), exhibits a dramatic changeover around e=1.3 GeV/fm 3 , 
saturating for higher e values, while it is falling below. The syst. error on q is 
estimated to be ~ 30%. It is assumed that the number of nucleons participating in 
the collision (N) is proportional to the volume of the particle source at the thermal 
freeze out ||. The new results from Si+Au at 14.6 A GeV and p+p at 158 A GeV 
shown in figure 1, which are not included in jg], have been estimated using data from 
and methods described in Furthermore, px rises with N respectively with V 
below e=1.3 GeV/fm 3 while it does not depend on N respectively on V above e=1.3 
GeV/fm 3 . To illustrate this, two values of V are noted on figure 1. The changes 
of and 7r ± with N within the Pb+Pb system, have been first realized in ]10| . A 
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Figure 1. Initial energy density (e) dependence of: (a) The K + multiplicity over 
the effective volume of the particle source at thermal freeze out. (b) The J/fy/DD 
(measured/'expected') ratio H. 
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Figure 2. The parameter a, resulting from the N a fit to hadron yields shown 
as a function of (a): the mass of the particles in the region e >1.3GeV/fm 3 at 
SPS and (b) of the ^/s, for kaons. N is the number of participating nucleons. 



similar behaviour as the one seen in figure 1, can be inferred for pions as well as for the 
K/ir ratio (S.K. work in progress). The N a exponent of hadrons with (u,d,s) quarks 
above e=1.3 GeV/fm 3 , do not depend on the particle mass (figure 2, (a)). At e >1.3 
GeV/fm 3 a is near to one, as expected in case of a chemically equilibrated state, 
assuming N ~ V. The deviations seen in <fi, tt° and p may be due to the transverse 
momentum acceptance. Therefore, figure 2 (a) supports the assumption of a high 
degree of chemical equilibrium reached above e=1.3 GeV/fm 3 , among hadrons with 
u,d,s quarks. The N a exponent of kaons is found to depend strongly on ^/s for kaons 
(figure 2,(b)). Therefore, below e=1.3 GeV/fm 3 , pk (figure 1 and figure 2 (b)), and 
the K/t: ratio, show an increase with increasing N respectively with V. 
Figures 1 and 2 can be interpreted in two ways. Firstly, kaons may achieve a higher 
degree of chemical equilibrium only for e > 1.3 GeV/fm 3 , and may not be fully 
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equilibrated below [gj. The equilibration of strangeness is expected in a QGP and 
its observation at e ~ 1.3 GeV/fm 3 could therefore be a sign of a transition to QGP. 
In this case, it is a transition from a non equilibrated hadron gas to an equilibrated 
QGP. It is therefore not a well defined phase transition in the thermodynamic sense. 
Secondly, kaons can be in fact chemically equilibrated also below e= 1.3 GeV/fm 3 , 
and the change respectively the constancy of pk with Vf and a observed in figure 
1, can be a result of the increase of the freeze out temperature with ei below e= 1.3 
GeV/fm 3 , respectively of the stability of T/ above 1.3 GeV/fm 3 . This dependence of 
Tf on 6i, namely rising until it reaches a critical T c value and saturating above for 
all reactions, would strongly support the QCD phase transition appearing at e ~ 1.3 
GeV/fm 3 . This interpretation fully agrees with thermal models which suggest that 
particle ratios at freeze out are compatible with thcrmalization even in A+A collisions 
at 1 A GeV || . However the first interpretation is not in gross disagreement with || , 
because there the thermal model description is modified (introducing e.g. pk ~ V) in 
order to describe the data at 1 A GeV. 

Furthermore, the correct interpretation can be corroborated by further investigations 
discussed in the following. The nonzero baryochemical potential (/zb), which in the 
reactions shown in figure 1, happens to change with ti, makes the intepretation of 
figure 1 difficult. Therefore, it appears that the dependence of the temperature at 
chemical freeze out extrapolated to /if,=0, on ej, would help to identify and prove the 
QCD phase transition. A rising and then a for ever saturating freeze out temperature 
above e 1.3 GeV/fm 3 is a strong argument that the QCD phase transition occurs at 
this e, and figure 1 is a direct consequence of it. 

The question if the QCD phase transition appears at the critical in any volume, or 
if there is additionally a critical initial volume of the particle source above which the 
transition takes place, can be answered comparing QGP signatures in systems with 
different volumes but the same e^. For example comparing p + p, e + e~ etc collisions 
to heavy ion collisions e.g. at the same e. This is not yet done for the signature of 
the J/\P suppression and it has to be clarified e.g. using Tevatron data ||. For the 
signature of strangeness enhancement it is suggested by figure 1 in || that there is 
indeed a critical initial volume, only above which strangeness is enhanced over p+p 
at the same q. This conclusion follows, if we assume that Tevatron reaches at least 
€i values similar to SPS A+A collisions jy} and if figure 1 in || is not biased by the 
model calculation 

If strangeness is indeed not equilibrated at e < 1.3 GeV/fm 3 , this may explain the 
decrease of the double ratio (-?T/7r)(A+A/p+p) with increasing y/s. In particular, a 
larger strangeness annihilation is enforced by equilibrium at SPS reducing the strange 
particle yield. However the assumption of non equilibrium of ss at low e is not 
nessecary here, since the above observation can be possibly traced back to e.g. the 
variation of ps{A + A)/p,s{p + p) with y/s in A+A collisions. Furthermore, in the 
context of QGP formation, it seems irrelevant to discuss e.g. ss enhancement in A+B 
over p+p collisions in a nonequilibrium situation. It is the very establishment of 
equilibrium in the (u,d,s) sector, which can reveal informations on QGP. 
The kaon number densities in p+p and A+B collisions in figure 1, (a) are similar, 
when compared at the same e^. See also [ fl2| for a discussion of universality of pion 
phase space densities. 

Our prediction for the N dependence of hadrons at RHIC and LHC is the N 1 thermal 
limit, as long as hadron yields are dominated by low transverse momentum particles. 
Furthermore, if the changeover of pk at e =1.3 GeV/fm 3 shown in figure 1 is due to 



Hadronic centrality dependence in nuclear collisions 



4 



the QCD phase transition, we predict for RHIC and LHC the same total strangeness 
(or kaon) number density and the same freeze out temperature, -after correction for 
the hb dependence-, as for e =1.3-3.0 GeV/fm 3 . If this change is however due to 
the onset of equilibrium in a hadronic gas, and the QCD phase transition takes place 
at higher e, it may manifest itself through a second changeover of hadron number 
densities, ratios and freeze out temperatures -after correction for the different \ib- e.g. 
in RHIC above e - 3 GeV/fm 3 . 

Assuming that the IMR dimuon enhancement seen by NA50 is due to open charm, 
the following observations can be made: a) open charm appears not to be equilibrated 
(a = 1.7) (figure 2, (a)) ||. b) The J/f/DD ratio deviates from p+p and p+A 
data also in S+U collisions (figure 1, (b)), above e ~1 GeV/fm 3 . c) It therefore 
appears that both charm and strangeness show a discontinuity near the same e~l 
GeV/fm 3 §, similar to the critical e c ~l-2 GeV/fm 3 predicted by QCD § §. d) 
The N dependence of the Jl^ / DD ratio can be interpreted as the J/ 1 ^ being formed 
through c,c coalescence e) Finally, the enhancement factors of hadrons with 
u,d,s,c quarks may be connected in a simple way to the mass gain of these particles 
in the quark gluon plasma (table below) T q are the enhancement factors of the 
lightest mesons with u,d,s,c quarks (tt, K, D), if they are produced out of a quark gluon 
plasma (e.g. g + g — > s + s (1)), as compared to their direct production from hadron 
interactions away from the transition point (e.g. p + p — > K + + A + p (2)). The gain 
is taken proportional to m part i C i e — m quar k s , as this expresses the different thresholds 
of reactions (1) and (2). In the table below the predicted enhancement factors (T q ) of 
hadrons with u,d,s,c quarks from a QGP are compared to the experimentally measured 
ones (Eg), and are found to be similar. (Definitions: th q = mo — m q , m u .d—7 MeV, 
m s =175 MeV, m c =1.25 GeV, m = m{ir, K, D)). 



Quark flavour th q T q = ^th q /th u 2 E ( N + N) E q = E/E Uid 

u,d 133 1 tt/N ~ 1.12 1 

s 320 1.55 K/N ~ 2 1.79 

c 615 2.15 DD meas/exp - 3 2.68 
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